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Automated tuning of bifurcations via feedback 
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1 Introduction 



u. 

' The present paper studies a feedback regulation problem that arises in at least 

<-j ■ two different biological applications. The feedback regulation problem under 

consideration may be interpreted as an adaptive control problem, but has not 
yet been studied in the control literature. The goal of the paper is to introduce 
this problem and to present some first results. 

The feedback regulation problem is concerned with a forced dynamical sys- 
tem 

> 

^ ■ i = U(x,u(t)) (1) 



which depends on a parameter fj,. The input u(t) in equation (Q) represents an 
' external stimulus; it is not a control variable. We are interested in the input- 

£SJ , to-state properties of equation (||). In general, the input-to-state properties 

' of equation (Q) depend on the parameter /i. In particular, the input-to-state 

properties may change drastically as a function of /i when the unforced dynamics 
exhibit a bifurcation. In this case, quite interesting and remarkable amplification 

S properties may result from operating in the proximity of the bifurcation point . 
We illustrate this with two biological examples from the literature. 
, The first example is concerned with the auditory system. In order to detect 

the sounds of the outside world, hair cells in the cochlea operate as nanosen- 
sors which transform acoustic stimuli into electric signals. In a series of recent 
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papers 0, ||, |[ pj, the hair cells in the cochlea are modeled as active, almost 
self-oscillating systems. Ions such as Ca ++ are believed to contribute to the 
hair cell's tendency to exhibit spontaneous oscillations. For low concentrations 
of the ions, the viscous damping forces of the fluid that surrounds the hair cells 
dominate and the hair cell oscillations are damped. As the concentration in- 
creases the system undergoes a Hopf bifurcation, the dynamics become unstable 
and the hair cells exhibit spontaneous oscillations. It is argued in these papers 
that the hair cells operate in the proximity of this bifurcation point, where 
the activity of the ions compensates the damping effects. In this case, even a 
weak sound stimulus can cause a detectable oscillation. (This follows from the 
generic properties of a forced dynamical system exhibiting a Hopf bifurcation) . 
This mechanism thus provides an explanation for the experimentally observed 
ultrahigh sensitivity of the ear. 

The second example arises in the study of persistent neural activity |j, |l2], 
H, H). Neural activity of a single neuron has a natural tendency to decay with 
a relaxation time of about 5-100 ms. This natural tendency to decay can be 
opposed by positive synaptic feedback loops. If this feedback is weak, then the 
natural tendency to decay dominates and neural activity still decreases. As the 
feedback gain is increased, the neural dynamics undergo a bifurcation and the 
dynamics become unstable. When the feedback is tuned to exactly balance the 
decay, then neural activity neither increases nor decreases but persists without 
change. According to a long-standing hypothesis, this is the mechanism that so- 
called neural integrators use to maintain persistent neural activity. A transient 
stimulus of a neural integrator can then cause a persistent change in neural 
activity. This mechanism forms the basis for short-term analogue memory and 
plays a central role in the oculomotor control system. 

Both examples illustrate the remarkable and interesting input-to-state prop- 
erties of a dynamical system poised at a bifurcation. It is clear, however, that 
operating in the proximity of a bifurcation point requires a fine-tuning of param- 
eters. And thus the question arises as to how a (biological) system can be tuned 
with high precision to its critical bifurcation point. In both examples it has been 
suggested that feedback regulation of the bifurcation parameter may provide a 
robust mechanism to ensure the required fine-tuning of parameters Q [|, [|, [l6| . 
In the literature on hearing this has given rise to the terminology of "self-tuned 
Hopf bifurcation" . Quite remarkably both studies have been pursued indepen- 
dently of each other in spite of their strong similarities which are obvious from 
a control engineering perspective. 

The present paper aims to initiate a mathematical study of this feedback 
regulation problem. The problem may be formulated as follows. Find an adap- 
tation law which steers the bifurcation parameter \x to its critical value /io 
without prior knowledge of this critical value /x . Our main motivation for 
studying this problem stems from the fact that the precise biophysical nature 
of the feedback mechanisms involved in the above biological applications is still 
unknown. A theoretical study of this problem may guide the search for pos- 
sible biophysical mechanisms. Even when the precise biophysical adaptation 
mechanisms will have been discovered, it is to be expected that a profound un- 



2 



derstanding of the mechanisms involved can only be achieved when insight has 
been gained in the fundamental trade-offs and limitations that may be inherent 
to this problem. This is very similar in spirit to the internal model principle 
which has turned out to be central to the understanding of adaptation in bacte- 
rial chemotaxis jl7j , and it may be seen as an illustration of the important role 
that control engineering may play in the emerging field of systems biology. 

A second motivation for the present study is of a more fundamental nature. 
In view of the quite remarkable and interesting input-to-state properties that 
arise when operating in the proximity of a bifurcation point, the problem under 
consideration may turn out to be of independent interest. Despite its mathemat- 
ical appeal, it has not attracted attention before in the control community. The 
present paper brings together several well-established techniques from nonlinear 
and adaptive control to derive some first results for this problem. 

We end this introduction with some references to related work. The present 
research bears some similarity with the problem of experimental instability de- 
tection where an operating parameter is adapted online in order to effi- 
ciently locate bifurcations through experiments. A different problem which is 
related to the present study concerns the detection and prediction of instabili- 
ties via closed-loop monitoring techniques || [To| . The possible implications of 
the present work for the problems studied in these papers remains a topic for 
further research. 



2 Self-tuning of a first-order system 

The one-dimensional system 

x = (jj, - no)x + u(t) (2) 

captures some of the essential features of the neural integrator. In this interpre- 
tation fiQ represents the natural decay rate of neural activity and /i corresponds 
to the synaptic feedback gain. We view /i as an adjustable parameter and /io as 
an unknown constant. Clearly this system exhibits a bifurcation. If /i = /io then 
equation @ behaves as a perfect integrator, if fj, < [iq (respectively /i > /io) 
then equation (||) is referred to as a leaky (respectively unstable) integrator. We 
ask the question as to how proximity to the bifurcation point may be ensured. 

We study this question from an engineering perspective. In a first approxi- 
mation we ignore the presence of the input and ask the following simpler ques- 
tion. Find an adaptation law for the parameter /i which steers /i to its bifurca- 
tion value /io for the system 

x = (n - fi a )x. (3) 

This adaptation law may depend on x and /i, but should be independent of /io, 
as this value is not known (precisely). Let us first discuss the feasibility of this 
problem. It is easy to see that if x = at some time instant then x — for 
all times and in this case it is clearly impossible to steer /i to /io without prior 
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knowledge of fio- We therefore restrict attention to the set of strictly positive 
values for x, which is invariant under the dynamics (^). Considering only strictly 
positive values for x is physically relevant, as this variable represents a level of 
neural activity (rate of action potential firing) . 

The following theorem provides sufficient conditions for the adaptation law 

A = fix) - g{n) (4) 

to steer fi to its bifurcation value fio for the system ([|). 

Theorem 1. Let (1q G R and consider continuously differentiable functions f : 
R>o — > R and g : R — > R. ylsswme i/iai / is strictly decreasing, g is strictly 
increasing, and g(/J,o) is in the image of f. Then the nonlinear system 
(g) with x € R>o /ias a unique equilibrium point, which is globally 

asymptotically stable (and locally exponentially stable ifdf/dx takes only strictly 
negative values and dg/d/j, only strictly positive values). At equilibrium /i equals 
Mo- 

The proof relies on a coordinate transformation which converts (||)-(g) into 
a nonlinear mass-spring-damper system. Global asymptotic stability follows 
readily from LaSalle's theorem. Local exponential stability is shown by means 
of the linearization principle. 

Proof. We introduce new coordinates q = ln(x) — ln(/ _1 ((;(/io))) andp = fx—fio- 
This transformation from (x,fx) to (q,p) is a global C°°-diffeomorphism from 
R>o x R to R 2 . Expressed in these new coordinates becomes 

q=P, (5) 
P = /(exp(g)/ _1 (.9(Mo))) -g(p + IM))- (6) 

The system of equations (||)-(||) has a unique equilibrium point at the origin. 
This equilibrium corresponds to an equilibrium in original coordinates where 
(j, = hq. 

First we prove that the null-solution of (|5j)-@ is globally asymptotically 
stable. We rewrite (||)-(§) as 

q = P, (7) 
P=-f(q)-9(p), (8) 

with /(g) = -f(exp(q)f- 1 (g(^o)))+g(fi ) and g(p) = g(p+^ ) -g{fMy)- Clearly 
/ and g are both strictly increasing, continuously differentiable functions from 
R to R which are zero at zero. The candidate Lyapunov function 

V:R 3 ^R:(q,p)~ J* f(£)d£ + V - (9) 
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is twice continuously differentiable, positive definite and radially unbounded. 
Its time derivative along the solutions of (0)-(||) satisfies 



n*P) = -P§<r)\ <0 a * Plt0 ' do) 

1 = lip = o. 

Since the null-solution is the only solution of (R)-(||) along which p vanishes 
identically, global asymptotic stability of the origin follows from LaSalle's the- 
orem [|[ Corollary 3.2]. 

Next we prove that the equilibrium of (||)-(||) is locally exponentially stable 
if d//da; takes only strictly negative values and dg/dp, only strictly positive 
values. Since the transformation from (x, p) to (q,p) is a C°°-diffeomorphism it 
suffices to prove that the null-solution of (@)-@ is locally exponentially stable. 
This follows readily from the linearization principle. The linearization of (||)-(^|) 
around the origin is given by 

q = P, (11) 

v = ^(/^(.g^o)))/- 1 ^^))? - (12) 

If d//dx takes only strictly negative values and dg/dp only strictly positive 
values then the linearization (|Il|)-(|l2|) is exponentially stable, as required. □ 

Discussion of Theorem [I] 

1. A mathematically appealing adaptation law which satisfies the assump- 
tions of Theorem [l] for all possible values of po is given by 

p = — ahi(x) — bp (13) 

with a, b > 0. This adaptation law ensures convergence to the bifurcation 
point for any possible value of /Jo- It has the interesting property that the 
coordinate transformation introduced in the proof of Theorem [l] converts 
(H]) , ([II]) into the linear system 

q = p, (14) 
p = —aq — bp. (15) 

From a biological point of view however this adaptation law is probably 



not very relevant, among others because the right hand side of (13) is not 
bounded. A bounded adaptation law which could be more relevant for 
biological applications is for example 

A=l/(l + a ; 2 )-l/(l + cxp(- M )). (16) 

This bounded adaptation law also ensures convergence to the bifurcation 
point for any possible value of fj,o ■ 
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2. There is a subtle relationship between self-tuning of bifurcations and ideas 
from robust control theory. This relation is made explicit by equation (|J) , 
which represents an integrator and corresponds to integral action stud- 
ied in robust control theory. Although perhaps surprising at first sight, 
this relation may be understood by regarding the constant /io as an un- 
known perturbation acting on the system. It is well-known from robust 
control theory that (under appropriate conditions) rejecting or tracking 
an unknown, constant disturbance requires integral action. (This is a spe- 
cial case of the internal model principle — see JIt], |l5) for a discussion of 
the internal model principle from a systems biology perspective). The 
present manifestation of integral action, however, differs from traditional 
robust control situations because the integral action is not generated by 
the adaptation law but is inherent to the dynamics of system (H) itself. 

3. If an additional perturbation term sp(x,n,t) is considered 

x = (ju - Ho)x + ep{x,n,t), (17) 

we may conclude that the adaptation law (^) steers /i approximately to 
/io under the conditions of Theorem ^ provided the perturbation is small 
enough. 

Theorem 2. Let /io € K and consider continuously differentiable func- 
tions f : M>o — > K and g : R — > M. Assume that f is strictly decreas- 
ing, g is strictly increasing, and g(fJ-o) is in the image of f. Consider a 
strictly positive parameter e and a continuous function p : R>o X R 2 — > 
R : (x, /i, t) — ► p(x, /!, t) and assume that p is bounded in t uniformly with 
respect to (x, /i) belonging to compact subsets o/M>o X M. Then the unique 
equilibrium point for the original system (^)^(^) is uniformly semiglobally 
practically asymptotically stable for the perturbed system ( |l7j ) -(^) . 

(We refer to the Appendix for a definition of uniform semiglobal practical 
asymptotic stability.) In particular, this theorem guarantees that along 
all trajectories of (fj~7])-([|) originating in an arbitrarily large compact sub- 
set of R>o x R, /i will eventually converge to an arbitrarily small open 
neighborhood of /io provided e is small enough. This result is a manifesta- 
tion of the well-known robustness of asymptotic stability with respect to 
small perturbations. It provides some justification for ignoring the exter- 
nal stimulus in equation (^) and studying the simpler equation (||) instead, 
at least in a first approximation. 



3 A second-order system 

In a series of articles ^, ^] it is argued that various nonlinear phenomena 
in the auditory system (such as ultrahigh sensitivity to weak signals) may be 
understood in terms of the generic properties of a forced dynamical system 
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exhibiting a Hopf bifurcation. In H this is illustrated by means of a standard 
model for nonlinear oscillations which (with the notation of the present paper) 
takes the form 

x + Oo - v)x + Ai 3 + u?x = u(t). (18) 

As before \i is an adjustable parameter and /io is an unknown constant. In the 
absence of an external input u, equation ( |l8| ) describes damped oscillations if 
(i < (Xq. At fi — fi Q the system undergoes a Hopf bifurcation and for /i > /i 
spontaneous oscillations are being generated. Assuming that it is possible to 
tune (j, = no such a model captures the biophysical nature of hair cell oscillations 
within the cochlea where the hair cells are thought of as active, almost self- 
oscillating systems rather than passive oscillators. In ||, |[ |tJ it is shown how 
such a dynamical system operating at a Hopf bifurcation indeed may account 
for several of the observed nonlinear phenomena in hearing. Again the question 
arises as to how proximity to the bifurcation point may be ensured. In [^j, [| it 
is suggested that a feedback mechanism is responsible for this self-tuning and 
numerical simulations are provided to support this claim. 

Here we want to contribute to a mathematical analysis of self-tuning of an 
oscillator. In a first approximation we ignore the cubic damping term and the 
external input and we study the following question. Find an adaptation law 
for the parameter \x which does not depend on fio and which steers /i to its 
bifurcation value (1q for the system 

x + (/Jo — fj)x + lo 2 x = 0. (19) 

As before let us first discuss the feasibility of this problem. It is easy to see 
that if x = x = at some time instant then x = x = for all times and in this 
case it is clearly impossible to steer /j, to /io without prior knowledge of (j,q. We 
therefore restrict attention to the set where x and x are not both zero — this set 
is invariant under the dynamics (|l9|). 

The following theorem asserts that the adaptation law 

fx = -a \n(y/x 2 + {x/uj) 2 ) - bfi (20) 

steers /i to its bifurcation value /j,q for the system (|l9|) if < a < b 2 and b > 0. 

Theorem 3. Let 6 R anda,b,uj G R>o- Assume that a < b 2 . The nonlinear 
system (fl9|) -(|2C|) with (x,x) G K 2 \ {(0,0)} and has a unique periodic 

orbi^\ which is globally asymptotically stable and locally exponentially stable. On 
this periodic orbit fi = /io. 

The proof of global asymptotic stability is based on passivity techniques. 
We interpret ©-© ( in a different coordinate system) as a feedback inter- 
connection of a linear system with a dynamic feedback which satisfies a sector 
condition. Via the Kalman-Yakubovich-Popov lemma we obtain a proper, non- 
increasing Lyapunov function. Global asymptotic stability of the periodic orbit 
then follows from LaSalle's invariance principle. 

A periodic orbit is a subset of the state-space which is the image of a periodic solution. 
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Proof. First we introduce new coordinates r and <j> according to x = rcos(0) 
and x = — ra;sin(^>). The transformation from (x,x) to (r, (f>) is a global 
C°°-diffeomorphism from R 2 \ {(0,0)} to R>o x 5 1 . Expressed in the coor- 



dinates (r, </>, /x), equations (19)~p0|) become 

r = (p - /x )r sin 2 (</>), (21) 

= w + (/i - ^o) sin(0) cos(</>), (22) 

/x = — aln(r) — bfi. (23) 

Next we introduce new coordinates q = ln(r) + 6/xo/a and p = /x — /xo- The 
transformation from (r, /x) to (g,p) is a global C°°-diffeomorphism from R>o x R 
to W 2 . Expressed in the coordinates (q,<fi,p), equations ([H|)-(|23|) become 



<? = psin 2 (0), (24) 
<fi = uj + p sin((f)) cos(cf)) , (25) 
p = —aq — bp. (26) 

Clearly the system of equations (p4])-(p6"|) has a periodic orbit {(q,<f>,p) G R x 
S 1 x R : q = p = 0} which we denote by A. This periodic orbit corresponds 
to a periodic orbit in original coordinates where /x = [1q. It is clear from the 
following paragraph that this periodic orbit is unique. 

First we prove that the periodic orbit A of (|24])-(|26]) is globally asymptot- 
ically stable. We interpret (|2^)-(|2^) as a feedback interconnection of a linear 
control system 

q = -u, (27) 
p=-aq- bp, (28) 

with a negative, dynamic feedback u = — sin 2 (c/))p where 4> satisfies (p5|). The 



transfer function H(s) of (27)— (128) from u to p is given by 



= s(sTbj' (29) 

Since a and b are strictly positive and a < b 2 the transfer function H(s) + 1 is 
positive real. Hence by (a modified version of) the Kalman-Yakubovich-Popov 
lemma (proven in and described in d, Exercise 10.2]) there exists a positive 
definite, quadratic Lyapunov function {q,p) > V{q,p) whose time derivative 
along the solutions of the control system (|27|)-(f28|) satisfies 

V(q,p,u) < up + u 2 . (30) 



In other words the control system (|27|)-(|28|) with input u and output p is input 
feedforward passive with a shortage of passivity jll| . Since the feedback satisfies 

up + u 2 = (sin 4 (0) - sin 2 (0))p 2 = -(sin(<£) cos(0)p) 2 , (31) 
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it follows that the time-derivative of V along the solutions of ([24|)-(|26|) satisfies 

V(q,<l>,p)< . (32) 

1 = it sm(0J cos(0)p = 0. 

Since (q, <f),p) t— > V(q,p) is positive definite with respect to the periodic orbit A 
it follows that A is stable. In addition, since (q,(j),p) i— > U(g,p) is radially 
unbounded with respect to A and since A is compact it follows that every 
solution of j2^)-(|26|) is bounded. Finally, since A is the largest invariant set of 
©-© contained in {{q,(f),p) elxS'xl: V(q,<j>,p) = 0} it follows from 
LaSalle's theorem that every solution of (p4|)-(p6|) converges to A. 

It remains to be proven that the periodic orbit of ([L9])-(|2(]) is locally ex- 
ponentially stable. This follows directly from Theorem |I|which is stated and 
proven below. □ 

Although mathematically appealing, the adaptation law (|2^) is probably not 
very relevant from a biological point of view, among others because the right 
hand side is not bounded. It would therefore be interesting to have a result 
available that applies to more general adaptation laws 



A = /(V^ + (i/w) 2 )-flO*). (33) 

This is the subject of Theorem |[ Unlike the previous results, Theorem ||] is a 
local stability result. 

Theorem 4. Let no 6l and consider continuously differentiable functions f : 
M>o — * M and g : M — > IR. Assume that g((J>o) is in the image of f. Consider 
r* e f^ 1 {g{^o)) and assume that < -(d//dr) (r*)r* < ((cWd^X^o)) 2 and 
(dg/d/i)(/i ) > 0. Then the system of equations (|l9|) and ( |33| ) with (x,x) S 
K 2 \ {(0,0)} and has a periodic orbit which is locally exponentially stable 

and where /i = /i - 

Proof. We proceed along the lines of the proof of Theorem ||. First we introduce 
new coordinates r and <f> according to x = rcos(<f>) and x = — rtv sm((f>) . The 
transformation from (x, x) to (r, (f>) is a global C°°-diffeomorphism from M 2 \ 
{ (0 , 0)} to R>o x S 1 . Expressed in the coordinates (r, <fi, /«), equations ( p^| ) and 
( P3[ ) become 

r = (fi- no)r sin 2 (0), (34) 

4> = oj + (p - ^o) sm(4>) cos((f>), (35) 

A = f(r) - g(ji). (36) 

Next we introduce new coordinates q = ln(r) — ln(r*) and p = ji — fio- The 
transformation from (r, fj,) to ((?,p) is a global C°°-diffco:rnorphism from M>o x K 
to M 2 . Expressed in the coordinates (q,<j),p), equations (^4|)-(|36l) become 

q = psin 2 {<!>), (37) 

= uj + psin(<j>) cos((/)), (38) 

P= /(exp(g)r*) -flf(p + Mo). (39) 
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Clearly the system of equations (p7|)~(p9|) has a periodic orbit {(q,4>,p) elx 
S 1 x R : q = p = 0} which we denote by A. This periodic orbit corresponds to 
a periodic orbit in original coordinates where fi — 

Since the transformation from (x, i, //) to (9, </>,p) is a C^-diffeomorphism it 
suffices to prove that the periodic orbit A of (|37|)-(|39|) is locally exponentially 
stable. We first ignore the "higher order" terms in the right hand side of 
and ( |39| ) and consider the simpler system 

q=psm 2 (<j)), (40) 
4> = u, (41) 
p = -aq - bp, (42) 

where we have introduced the notation a = —(df/dr) (r*)r* and b = (dg/d^)(/uo). 
By the assumptions of the theorem < a < b 2 and b > 0. Repeating the ar- 
guments of the proof of Theorem ^ it is easy to see that the set A is a globally 
asymptotically stable periodic orbit of ([i"o|)-([42"|). Because of the special struc- 
ture of equation (^1) we may interpret ^ as a time variable and the system of 
equations ( ^0| ) and ( 4^) as a periodically time- varying linear system. With this 



interpretation the periodic orbit A of (40)- (H2j) corresponds to the null-solution 



of the linear system of equations ( flop and (|42|). Standard converse theorems 
for periodically time-varying linear systems (see for example j^, Theorem 3.10]) 
yield the existence of strictly positive real numbers c\ , C2 and C3 and a continu- 
ously differentiable Lyapunov function (q,tfi,p) 1— > V(q, <fi,p) which is quadratic 
in (q,p) and (together with its partial derivative dV/d<p) bounded in <fi such 
that 

ci(q 2 +P 2 ) < V{q,cj>,p) < c 2 (q 2 +p 2 ) (43) 

and such that the time derivative of V along the solutions of (|4(]) (^) satisfies 

V(q,4>,p) < -c 3 (q 2 +p 2 ). (44) 

With this Lyapunov function V we now prove that the periodic orbit A of (|37|)- 
( ^9| ) is locally exponentially stable. Indeed, for q and p sufficiently close to zero 
the time derivative of V evaluated along the solutions of (|37|)-(j39l) satisfies 

V{qA,p) <-j(q 2 +p 2 ) (45) 

since the extra terms in the right hand side of J38| ) and ( |39| ) give rise to extra 
terms in the Lyapunov balance which are bounded in <f> and of order higher 
than two in (q,p). This shows that the periodic orbit A of (|37|)-(|3"9|) is locally 
exponentially stable. □ 

Discussion of Theorems |3| and [i] 

1. The assumptions of Theorem ^ involve the unknown critical value no- 
If /io is known to belong to some interval, it is of interest to have an 
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adaptation law which satisfies the assumptions of the theorem for all /io 
in this interval. For example, it is easily verified that the adaptation law 



1/(1 + (Vx 2 + {x/lo) 2 ) °) - bn (46) 



with < a < 4b 2 and b > satisfies all assumptions of the theorem if 
< fi < 



2. The proposed adaptation laws (|20|) and (|33|) depend not only on x and 
fi but also on x and to. Their implementation requires that the state 
variable x is measured and the parameter lo is known. This may be a 
disadvantage. 

3. If an additional perturbation term ep{x, x, fi, t) is considered 

x + (/i — /i)£ + u> x = ep(x, x, fj,, t), (47) 



we may conclude that the adaptation law ( j20| ) or (33) steers (j, approxi- 
mately to /io under the conditions of Theorem |^ or Theorem ^ provided 
the perturbation is small enough. 

Theorem 5. Let /io £ K a^d consider continuously differentiable func- 
tions f : R>o — > K a«c? g : R — > R. ylss«me that (/(/io) is m £/ie image 
of f . Consider r* 6 /~ 1 (ff(/io)) a^rf assume that < — (d//dr)(r*)r* < 
((d<?/d/i)(/xo)) 2 and (dg/d/i)(/io) > 0. Consider a strictly positive param- 
eter e and a continuous function p : (R 2 \ {0}) x R 2 — > R : (x,i,/i,t) — > 
p(x,x, n,t) and assume that p is bounded in t uniformly with respect 
to (x,x,fi) belonging to compact subsets of (R 2 \ {0}) x R. Then there 
is a subset of (R 2 \ {0}) x R where /i = /io and which is uniformly practi- 
cally asymptotically stable for the perturbed system j4^ ) and (|33|). If /(■) 
and c/(/z) iafce f/ie particular form — aln(-) respectively bfj,, then there is a 
subset of (R 2 \ {0}) x R w/sere fj, = )1q and which is uniformly semiglobally 
practically asymptotically stable. 

(We refer to the Appendix for a definition of the notion of uniform (semiglobal) 
practical asymptotic stability.) This result is a manifestation of the well- 
known robustness of asymptotic (or exponential) stability with respect to 
small perturbations. It provides some justification for ignoring the cu- 
bic nonlinearity and external stimulus in equation (^|) and studying the 
simpler equation (J19) instead, at least in a first approximation. 
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A Practical and semiglobal stability definitions 



Let fl be an open subset of R™ (n £ N) and consider a family of continuous 
functions f £ : K x SI — > R™, labeled by a parameter e > 0. We are interested in 
the stability properties of the family of differential equations 

x = f s (t,x) (48) 

for small values of e. Let A be a compact subset of f2, which need not necessarily 
be forward invariant for (|48|). 

Definition 1. For the family of differential equations ([48|) ; the set A is: 

1. Uniformly practically stable if for every open neighborhood XJi of A there 
is e* > and an open neighborhood U\ of A such that for all e £ (0, £*] 
every solution £ of ( f48| ) satisfies: if £(io) € £A /or some to w the domain 
of then £(t) £ U2 for all t >to in the domain of £. 

2. Uniformly practically asymptotically stable if it is uniformly practically 
stable and, in addition, there is an open neighborhood U\ of A such that 
for every open neighborhood U2 of A there is T > and e* > such 
that for all e £ (0, e*] every solution £ of ( f48| ) satisfies: if £(io) £ C^i /° r 
some to in the domain of then £(t) £ U2 for allt > to + T in the domain 
oft 

3. Uniformly semiglobally bounded if for every compact subset K\ of CI there 
is e* > and a compact subset K2 of fl such that for all e £ (0, e*] every 



solution £ of (48) satisfies: if £(to) £ K\ for some to in the domain of t;, 



then £ K2 for all t >to in the domain of £. 

Uniformly semiglobally practically asymptotically stable if it is uniformly 
practically stable and uniformly semiglobally bounded and, in addition, for 
every compact subset K of CI and every open neighborhood U of A there 
is T > and e* > such that for all £ £ (0, £*] every solution £ of ( [48| ) 
satisfies: i/£(io) £ K for some to in the domain of t;, then £(t) £ U for 
all t > to + T in the domain of £. 
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